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Abstract 

It is shown that, for isolated many-electron Coulomb systems with Coulombic 
external potentials, the usual reductio ad absurdum proof of the Hohenberg-Kohn 
theorem is unsatisfactory since the to-be-refuted assumption made about the 
one-electron densities and the assumption about the external potentials are not 
compatible with the Kato cusp condition. The theorem is, however, provable 
by more sophisticated means and it is shown here that the Kato cusp condition 
actually leads to a satisfactory proof. 
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1. Introduction 

The Hohenberg-Kohn theorem [1] underlies the foundation of the density functional 
theory [2] and since 1964 when it was formulated and proved by reductio ad absurdum, 
it has significantly influenced the state of art of quantum theory of atoms, molecules, 
clusters, and solids. 

The aim of the present work is an attempt to rethink and reanalyze the original 
proof by reductio ad absurdum of the Hohenberg-Kohn theorem. It is shown here that, 
although the result is a generally correct one, the original proof cannot be maintained 
when the external potential is of Coulomb form because it implies that the supposed 
wave functions must violate the Kato electron-nuclear cusp conditions. However, more 
sophisticated proofs can be adduced, in particular directly from the Kato theorem, 
which avoid this problem. 

2. General Part 

On p. B864 of their work [1], Hohenberg and Kohn state that they "... develop an 
exact formal variational principle for the ground-state energy, in which the density" 
p(r) (in a widely accepted notation) "is the variable function. Into this principle enters 
a universal functional" F[p(r)], "which applies to all electronic systems in their ground 
state no matter what the external potential is." Following Hohenberg and Kohn [1], 
let us consider "a collection of an arbitrary number of electrons, enclosed in a large 
box and moving under the influence of an external potential v(r)" where r e 3? 3 "and 
mutual Coulomb repulsion." The Hamiltonian H„ of a given iV-electron system casts 
as 

H v N = T e N + V e N e +V N (1) 
where is the kinetic energy operator of N electrons, Vj£ is the corresponding inter- 
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electronic Coulomb operator, and 



V 



N 




(2) 



is the total external potential. Hohenberg and Kohn [1] further assume (p. B865) that 
possesses the least bound-state (ground-state) wavefunction ^ ( r i, r 2, £ 
?i 1 (3? 3Ar ) (spins are omitted for simplicity) and the latter is nondegenerate within the 
total coordinate - spin representation. H l (^ 3N ) is the Sobolev space of iV-electron 
wavefunctions whose norms defined as (/ n£Li d 3 Ti(\ ^ | 2 + | I 2 ) 1 / 2 < oo. Let us 
then define the corresponding ground-state one-electron density [3] 



where = { Po e L 1 ^ 3 ) | p > 0, ^fp~ e H 1 ^ 3 ), J d 3 r Po (r) = N}, "which is clearly 
a functional of v(r) n (p. B865, Ref. [1]) since H„ is explicitly determined by v(r) 
under the fixed N,T^, and V^f. That is, there exist such mappings 



from L 3 / 2 (9ft 3 ) + L°°(!R 3 ) as the domain of external potentials [4, 5] (see precisely Eq. 
(2.4) in Ref. [6]) to 7i 1 ($R 3JV ) and further to if H„ does possess the nondegener- 
ate ground state. Notice that belonging of any v(r) to L 3 / 2 (K 3 ) + L°°(K 3 ) does not 
guarantee that the corresponding of the type (1) has this property [4]. Notice also 
that any Coulomb potential, v(r) = —Z a / | r — R a |, or any finite linear combination 
of Coulomb potentials, v(r) = -^ =1 Z a / | r - R Q |, belong to L 3 / 2 (K 3 ) + L°°(K 3 ) [5]. 
The mapping Cn in (4) is defined as that from L 3 / 2 (9? 3 ) + L°°(3? 3 ) as the subdomain 
of external potentials for which possesses the ground state into Ti l {$i 3N ). It is 
readily to prove that Cm is invertible (see, e. g., Ref. [5]) that is any pair of exter- 
nal potentials differed from each other by more than a constant determine a pair of 
different ground-state wavefunctions. The other mapping in (4) is the mapping T>n 




i=2 



(3) 



,(r)^ (r ll r 2 ,...,r,)^ /)o (r) 



(4) 
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from the ground-state iV-electron eigenwave functions onto the set C of the 
ground-state one-electron densities. An invertibility of the mapping V N relies on the 
Hohenberg-Kohn theorem (see also Eqs. (27)- (29) in Ref. [5]). 

Hohenberg-Kohn theorem [1]: u v(r) is a unique functional" of p(r), "apart 
from a trivial additive constant." 

Proof (p. B865, Ref. [1]): "The proof proceeds by reductio ad absurdum." We 
assume the existence of two "external" potentials i>i(r) and v 2 (r) such that 

t>i(r) 7^ v 2 (r) + constant. (5) 

Via Eqs. (2) and (1), t>i(r) and v 2 (r) define the Hamiltonians Hf and H 2 associated 
with two different iV-electron systems. Let us further assume the existence of the 
ground-state normalized wavefunctions *W G H l {^ N ) and ^ e H 1 ^ 31 *) of 
and H 2 , respectively. By virtue of Eq. (3), and yield the corresponding 
ground-state one-electron densities p£^(r) and p@\r). Hohenberg and Kohn [1] finally 
assume that 

(i ) tf(D ^ (2) 

(") pi 1) W=pi 2) (r) = Po(r). 

Applying the Rayleigh-Ritz variational principle, one obtains 

*(*)) 

V 2 N | ^} 

(6) 
(7) 
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= E( 2 )+ | ^r[ Vl (r)- V2 (r)] Po (r) 
= - j dh\ Vl (v) - v 2 (v)] Po (v) 



where the used formulas are indicated above the signs. 

Hohenberg and Kohn then conclude (p. B865, Ref. [1]) that adding (6) to (7) 
"leads to the inconsistency" 

eM+eW <EU+Ei 2 \ (8) 

and therefore, (8) implies that the assumption (ii) fails. "Thus v(r) is (to within a 
constant) a unique functional of p(r), "since, in turn, v(r) fixes" H„ "we see that the 
full many-particle ground state is unique functional of p(r). Q. E. D. 

Examine Eq. (8). It is obviously self-contradictory. (8) is deduced under the 
assumption that (5) is true together with the to-be-refuted assumptions (i) and (ii), 
both composing the negation of the statement of the Hohenberg-Kohn theorem. (8) 
then appears to be absurd in a sense of being obviously false and therefore the statement 
of the Hohenberg-Kohn theorem is correct. Logically speaking, the fact that Eq. (8) 
looks absurd implies that, first, one of the to-be-refuted assumptions, (i) or (ii), or 
simultaneously both, (i) and (ii), lead to the contradiction with (5) or, second, they 
are a priori invalid in a sense that one of them or both are somehow incompatible with 
(5). In the latter case, the statement of the Hohenberg-Kohn theorem is invalid unless 
it is proved in the other way. Note also that (i) and (ii) are obviously inconsistent for 
N — 1 and for the two-electron noninteracting systems. Explicitly, after deriving (8), 
one has to consider the following cases: 

(I) *W = *( 2 ) = * . 

This one directly yields pW = = p , that is, (ii) does hold. It also implies that 

(T N + V N )ty 

^ o 

if V-f and V<f are multiplicative operators, as suggested by Eq. (2). It is trivial to 
conclude that (9) contradicts to (5). However, there is no "inconsistency" because the 
last terms in the last lines of Eqs. (6) and (7) simply vanish. 
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(II) *(D^(2) ail d pW^fg). 

This case precisely lies in the line of the original arguments by Hohenberg and Kohn 
[1] proving thus that different external potentials determine different ground-state one- 
electron densities. 

(III) = *( 2 ) and pW±pf>. 

These two relations contradict to each other due to (3). 

(IV) A self-contradiction (ad absurdum) of Eq. (8) might also mean that the to- 
be-refuted assumptions (i) or/and (ii) in the original proof of the Hohenberg-Kohn 
theorem are self-contradictory with Eq. (5) and this is precisely the case of real many- 
electron Coulomb systems with Coulombic external potentials. In other words, the 
original reductio ad absurdum proof of the Hohenberg-Kohn theorem based on the 
assumption (5) is incompatible with the ad absurdum assumption (ii) due to the validity 
of the Kato theorem for such systems [6] . A similar statement is valid with regard to 
the proof by reductio ad absurdum of the invertibility of the mapping T>n in Ref. [5]. 

Let us recall that the Kato theorem [6] (see also Refs. [7, 3]) determines the 
character of the singularity of the exact A^-electron eigenwavefunction of at the 
electron-nucleus coalescences where the Coulombic external potential v(r) (see Eq. 
(2.2) and the conditions i) and ii) on p. 154 and Theorem I on p. 156 of Ref. [6]), 

"to = -^lu^W-r (10) 

(in atomic units) is singular. In Eq. (10), the ath nucleus with the nuclear charge Z a is 
placed at R a G 3? 3 . Any A^-electron eigenwave function \& of with v(r) of the form 
(10) and its one-electron density p$ then satisfy the electron-nucleus cusp conditions 

A^««( ri> r2) n_ 1: r i: r i+ i, rjv) \ n =n a = -Z a ^(r u r 2 , r^, R a , r m , t n ), 
an 

^p% v (r) \ r=Ra = -2Z a p*(R a ) (11) 
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where the superscript 'ai>' means the average of \I/ over an infinitesimally small sphere 
centered at and that of at r, and where i — 1, 2, N in the former relationship. 

Therefore, the true one-electron density of the given iV-electron system, moving 
in the Coulombic field of point nuclei, exhibits cusps (local maxima) at the positions 
of the nuclei. Assuming that the ground-state one-electron density p ( r ) is given and 
analyzing its topology over the whole coordinate space K 3 , one locates the positions of 
its cusps and evaluates there the lhs of the last equation (11). Altogether, the positions 
of the electron-nucleus cusps (as being always negative that follows from Eq. (11)) and 
the halves of the radial logarithmic derivatives of Pg V (r), taken with the opposite sign 
at these points, fully determine the Coulombic external potential v (r), Eq. (10), of the 
given system. That is, p ( r ) uniquely determines the external Coulombic potential v(r). 
Such interpretation of the Hohenberg-Kohn theorem for Coulombic external potentials 
was originally proposed by Coleman [8], Bamzai and Deb [9], Smith [10], and E. Bright 
Wilson (quoted by Lowdin [11]). Also, this naturally proves an invertibility of the 
mappings Cjv, V N , and V n Cn- Therefore, if a given pair of iV-electron systems with 
the Hamiltonians and of the type (1) are characterized by the same ground- 
state one-electron densities (that is equivalent to the to-be- refuted assumption (ii)), 
their external potentials vi(r) and v 2 (r) of the form (10) are identical, and thus and 
are identical as well. This contradicts (5) and hence, the assumption (ii) cannot 
be used, rigorously speaking, in the proof via reductio ad absurdum of the Hohenberg- 
Kohn theorem together with the assumption (5) for the Coulomb class of external 
potentials. In other words, (5) and (ii) are Kato-type incompatible to each other for 
purely Coulombic external potentials. 

Vice versa, the nuclei of the given A-electron system are isolated 3D point attractors 
behaving topologically as critical points of rank three and signature minus three [12]. 
However, there exist some "specific" many-electron systems whose ground-state one- 
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electron densities have local maxima at non-nuclear positions [13]. These local non- 
nuclear maxima might be the true ones or appear as a consequence of an incomplete, 
inadequate quantum mechanical treatment. Therefore, despite the present conclusion 
that in the original proof by reductio ad absurdum of the Hohenberg-Kohn theorem, 
one of its to-be-refuted assumption (ii) is incompatible, by virtue of the Kato theorem, 
with the assumption (5) 2 , the Kato theorem itself guarantees the existence of the one- 
to-one correspondence between the Coulomb class of external potentials (10) and the 
ground-state one-electron densities for nearly all many-electron except probably those 
"specific" ones. However, local non-nuclear maxima of the ground-state one-electron 
densities of these "specific" many-electron systems can be easily excluded since, first, 
they are not cusps at all and second, they disappear under some changes (i. e., bond 
lengthening) of the nuclear skeleton [13]. 

Summarizing, the Kato theorem implicitly contains the proof of the Hohenberg- 
Kohn theorem for iV-electron systems with the Coulomb class of external potentials 
and assures an invertibility of the one-to-one mapping V N C^ between the Coulombic 
external potentials and thus between many-electron Hamiltonian and the corre- 
sponding nondegenerate ground-state one-electron densities (see also Ref. [5]). Obvi- 
ously, due to its original formulation [6], the Kato theorem cannot be used for other 
types of external potentials (that is why the present treatment is only confined to 
the Coulombic ones of the type (10)) although it might be generalized to include the 
subclass of the singular ones, like, for example, a Yukawa potential that results in 
nonvanishing cusps (see p. 156 of Ref. [6] for the definition of the class of generalized 
Coulomb potentials). 

According to the work [1] of Hohenberg and Kohn, the Hohenberg-Kohn theorem 

implies the existence of the universal energy density functional for any isolated many- 

2 For a similar proof of the ensemble generalization of the Hohenberg-Kohn theorem see Section II 
of Ref. [14]. 
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electron Coulomb system. This statement has been usually interpreted as the second 
Hohenberg-Kohn theorem [2]. In the density functional theory, there exist two rigorous 
constructions of the universal energy density functionals based on their own rigorous 
proofs of the Hohenberg-Kohn theorem. This is the Levy-Lieb energy density func- 
tional [15, 4] and the energy density functional based on the group of the local-scaling 
deformations in 3? 3 consisting of topogical deformations mapping or topologically de- 
forming any pair of one-electron densities to each other [16, 2c]. The related Jacobian 
of such deformation gives rise to the first-order nonlinear differential equation, a so 
called "Jacobian equation" [17, 18] whose solution, namely, the corresponding defor- 
mation, does exist within this approach and it is unique [16, 2c]. Solving then the 
"Jacobian equation" enables to determine the deformation for any pair of well-behaved 
one-electron densities and to consistently extend the action of the local-scaling defor- 
mation group onto 7i 1 (K 3JV ) [16, 2c]. This larger local-scaling transformation group 
partitions 7i 1 (K 3JV ) into disjoint classes, orbits (see Refs. [16, 2c] for details). All orbits 
exhaust H 1 (?R. m ) and within an each given orbit, there establishes the one-to-one cor- 
respondence between its wavefunctions and the one-electron-densities. That is, these 
orbits are endowed with the characteristic that there are no two wavefunctions belong- 
ing to the same orbit that have the same density. Each orbit, say 0^ a \ is determined 
by its generator wavefunction vpj^l. Therefore, for a given orbit 0^ a \ one defines the 
energy density functional E a [p{r)] as merely the restriction of the energy functional 
E\$>\ on those wavefunctions that belong to the ath orbit, \P G O^. It is trivial to 
prove that, first, there are as many different energy density functionals as the orbits 
in 7i 1 (K 37V ), and second, that each density functional E a [p(r)] implicitly depends on 
the generator wavefunction vpj^. Evidently, the exact ground-state iV-electron eigen- 
wave function of the given Hamiltonian operator H belongs to a certain orbit, called 
the Hohenberg-Kohn one, 0^ HK ^ [16, 2c]. Within 0^ HK \ the Levy-Lieb energy den- 
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sity functional [15] exactly coincides with Ehk[p(y)] defined within the local-scaling 
deformation approach [16]. 

The explicit form of any energy density functional E a [p(r)} for any a has been ob- 
tained and the corresponding variational Euler-Lagrange equation has been also derived 
in Ref. [16] (see also Ch. 7 and 8 of Ref. [2c]). The rigorous mathematical framework 
of the local-scaling deformation approach to the density functional theory based on 
the "Jacobian equation" has recently been elaborated in Ref. [17]. The local-scaling 
deformation approach to the density functional theory has been also generalized for 
the spin densities, the momentum space representation, excited states, fractional occu- 
pation numbers, and finally to study nonadiabatic effects. The corresponding analogue 
of the Kohn-Sham approach has been also formulated in terms of the orbit generators 
(see Ch. 8 in Ref. [2c]). A variety of theoretical and computational applications of the 
local-scaling deformation density functional theory to atoms and molecules has recently 
been elaborated as well (see Ref. [19] and references therein). 
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